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^ • Abstract 
CD ! 

r~| , We analyze the contribution of the SUSY particles to the coupling of the lightest 

Higgs boson to two photons in supersymmetric theories. We discuss to what extent 
these contributions can be large enough to allow for a discrimination between the 
^ • lightest SUSY and the standard Higgs particles in the decoupling limit where all 

^ other Higgs bosons are very heavy and no supersymmetric particle has been dis- 

covered at future colliders. We find that only chargino and top squark loops can 
generate a sizeable difference between the standard and the SUSY Higgs-photon 
couplings. For masses above 250 GeV, the effect of chargino loops on the two- 
photon width is however smaller than ~ 10% in the entire SUSY parameter space. 
Top squarks heavier than 250 GeV can induce deviations larger than 10% only if 
their couplings to the Higgs boson are large. Since top squark contributions can 
be sizeable, we derive the two-loop QCD correction to squark loops and show that 
they are well under control. 
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1. Introduction 



In supersymmetric (SUSY) theories, the Higgs sector is extended to contain at least two 
isodoublets of scalar fields. In the minimal version, the Minimal Supersymmetric Standard 
Model (MSSM), this leads to the existence of five Higgs particles: two CP-even Higgs 
bosons h and if, a CP-odd or pseudoscalar Higgs boson A, and two charged Higgs particles 
Besides the four masses, two additional parameters are needed to describe the 
Higgs sector at tree-level: tan/3 the ratio of the two vacuum expectation values and 
a mixing angle a in the CP-even sector. However, only two of these parameters are 
independent, and choosing the pseudoscalar mass Ma and tan/3 as inputs, the structure 
of the MSSM Higgs sector is entirely determined. 

If the pseudoscalar mass Ma is very large, 3> M|, the pattern of Higgs masses is 
quite regular. The heavy CP-even, CP-odd and charged Higgs bosons are nearly mass 
degenerate, M^ ~ Mh± ~ M4, while the lightest CP-even h particle reaches its maximal 
mass value. At tree level, this value is simply a function of tan/3, M^^^ = Mz \ cos 2/3 1 < 
Mz- However when including the radiative corrections p|, which grow as the fourth 
power of the top mass and logarithmically with the common squark mass, the upper 
bound is shifted upwards, Mj^^^ ~ 130 GeV. In this so called decoupling limit [Q, which 
in practice is reached for Ma ~ 300 GeV, the lightest SUSY Higgs boson h has almost 
the same properties as the SM Higgs particle and the MSSM and SM Higgs sectors 
look practically the same, with one light Higgs boson with a mass below ~ 130 GeV. 

In the case where no genuine SUSY particle and no additional Higgs boson has been 
discovered at future high-energy coUiders, the task of discriminating between the lightest 
SUSY and the standard Higgs bosons, and therefore between the MSSM and the SM, in 
the decoupling limit is challenging. Indeed, since both have almost the same couplings to 
fermions and vector bosons, the production rates and the decay branching ratios [when 
SUSY Higgs decays are kinematically not allowed] are practically identical. 

Only indirectly that one can distinguish between the two models: if the SM is extended 
to the GUT scale, the value rrit — 175 GeV requires a Higgs boson heavier than Mfjo > 130 
GeV in order that the vacuum remains stable; since in the MSSM, M^ is constrained to 
be lighter than ~ 130 GeV, the measured Higgs mass will allow to discriminate between 
the SM and MSSM scenarios [^]. However, one can assume that new physics beyond the 
SM exists at a scale A < 10 TeV and in this case, Higgs masses in the range M^o ~ Mz 
will be still allowed. Furthermore, the SUSY Higgs mass bound M^ < 130 GeV is valid 
only in the MSSM: for more general SUSY scenarios where the Higgs sector is even more 
complicated [for instance in the NMSSM where an additional Higgs singlet is added], the 
upper bound on the lightest Higgs mass from triviality can be extended to M™'^^ ~ 150 
GeV [0, leaving a room for an overlap between the allowed h and masses. 
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A more "direct" way to discriminate between the standard and the hghtest SUSY 
Higgs particles is to look at loop induced Higgs boson couplings such as the (^gg [§], 
and $77 |jlO[ couplings, $ = or H^. In the SM, these couplings are mediated 



by heavy quark and W boson loops [only quark loops for the H^gg coupling]: since 
their couplings to the Higgs boson grow with the mass, they balance the decrease of the 
triangle amplitude with increasing loop mass, and these particles do not decouple even for 
masses much larger than Mho. In supersymmetric theories, additional contributions will 
be induced by loops with charged Higgs bosons, charginos and sfermions; Fig. 1. However, 
since the SUSY particles do not couple to the Higgs boson proportionally to their masses, 
their contributions are expected to be rather small for large masses. For very heavy SUSY 
particles, the loop induced vertices reduce to their SM values and again, no distinction 
between the SM and the MSSM can be made. 

The (^gg vertex can be measured in the main Higgs production process gg ^ ^ aX 
hadron colliders, or via the branching ratio BR($ gg) with the Higgs boson produced 
at e^e~ colliders. At the LHC, the determination of the cross section a{gg —>■ $) to 
the level of ten percent is rather difficult, due to uncertainties from the QCD corrections 
[which at next-to-leading order are very large |[Tl], [1^, increasing the cross section by 



almost a factor of two] and to a lesser extent from the parton densities. The branching 
ratio BR($ — > gg) is of the order of a few percent for M$ ~ 100 GeV, and its measurement 
at e~^e~ colliders with an accuracy of more than a few ten percent is also very difficult 
due to the contamination from charm and bottom quarks 



The $^7 vertex can be measured in the decay Z $7 at LEP and SLC if M$ < Mz, 
or in the reverse decay $ — >■ Z'y if M$ > Mz with the Higgs boson produced in the gg ^ ^ 
fusion mechanism at the LHC. However the rates are very small, BR(Z — > $7) < 10^^ 
and BR($ — > Zj — > < 10~^, leading to only a few events at LEP or the LHC and 

making the determination of the $^7 vertex with a reasonable accuracy very difficult. At 
future e+e~ colliders with the expected integrated luminosities / £ ~ 50 fb~^, running a 
few months on the Z resonance would allow to obtain a large sample of Z ^ $7 events if 



M$ < Mz] a precise measurement of the $^7 coupling would be possible in this case [^ 

The prospects for measuring the loop induced $77 vertex are as follows: 

[i) At the LHC the production rate for light Higgs bosons is very large, (T{gg ^ $) ~ 
100 pb ||15|, and despite of the small branching ratio BR($ 77) ~ 10~^, one would still 
have O{10^) 77 events after filtering out most of the background events, if the luminosity 
is high enough, C ~ 10^^ cm~^s"^. However, as discussed earlier, besides the uncertainties 
from the parton densities, the theoretical prediction of the production cross sections is 
affected by large uncertainties from higher QCD corrections. Since one measures only 
crx BR, a clean extraction of the $ — 77 width will be rather difficult. 
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(ii) At e~^e~ collider^, the main production mechanisms for Higgs particles are the 
bremsstrahlung process e~^e~ and the WW fusion process e"'"e~ W*W* $t'e^'e 

[p^ . At energies in the range of y/s ~ 500 GeV, the cross sections are around 100 fb for 
each process; even for integrated luminosities of ~ 50 fb~^ one would have only a few 
$ — > 77 events, a sample which does not allow a precise measurement. At higher energies 
the cross section for the WW fusion mechanism increases logarithmically: at ^/s ~ 1.5 
TeV and with J C ^ 200 fb~^, one would have (9(100) events allowing for a decent 
measurement. However, if no SUSY particles have been found at this energy, their effect 
on the $ — > 77 width will probably be too small to be visible. 

{Hi) The most promising way to have access to the $77 coupling is via the single 
Higgs production in the fusion process 77 ^ $ [18-22], with the photons generated by 



Compton-back scattering of laser light pSf . One can tune the energy of the 77 collider 
such as to produce the Higgs boson as a resonance in the s-channel. If the luminosity of 
the 77 collider is of the same order as the luminosity of the original e~^e~ collider, large 
production rates can be obtained. A measurement of the $ — > 77 partial decay width 
with a precision of the order of 10% could be feasible as will be discussed later. 



There are several studies of the Higgs-photon coupling in the MSSM |26[ which how- 
ever mainly focussed on the detectability of the /i — 77 signal at the LHC In this paper, 
we analyze this coupling with a different perspective: we scan the entire MSSM parameter 
space and single out the regions where the SUSY loops could give significant contribu- 
tions. Our aim is to answer to the important question [since the measurement of the 
Higgs-photon coupling is one of the most important goals of the presently discussed 77 
colliders] of how well one needs to measure the $ ^ 77 width in order to discriminate 
between the SM and the MSSM Higgs boson in the decoupling regime, if no SUSY par- 
ticle has been observed directly at the LHC or at an e~^e~ collider with a cm. energy of 

= 500 GeV. 

The paper is organized as follows. In the next section, we describe the Higgs sector in 
the decoupling limit and present for completeness the formulae for the loop contributions 
to the $77 couphng. In section 3, we analyze the Higgs production at 77 colliders, and 
estimate the precision with which the Higgs-photon coupling can be measured. In section 
4 we discuss the various contributions and isolate the parameter space in which these 
contributions are significant. Our conclusions will be given in section 5. In the Appendix, 
we derive the QCD correction to the squark loop contribution to the $77 amplitude. 



^At e^^e^ colliders one can also measure the e+e^ — > $7 cross section which is built up by loops of 
heavy particles [p^ ; however the cross sections are rather small, and large luminosities will be required. 
Another possibility is provided by the process 76" ^ e~$ as recently discussed in [0. 
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2. The Higgs-Photon coupling in the MSSM 

In the MSSM, using Ma and tan (3 as input parameters, and including the leading radiative 
correction which can be parameterized in terms of the quantity P] 



v^7r2 sin2/5 



(1) 



with TTiq the common squark mass, the CP-even Higgs boson masses are given by 



It 



\ 



1-4 



MlMl cos^ 213 + e{Ml sin^ /3 + M| cos^ /3) 



(2) 



In the decoupling limit, M\ ^ M|, the Higgs masses approach the values 



H 



M|cos2 2/5 + esin2/3 



1 + 



eMl cos^ (3 



Mfsin^ 2/3 + ecos2/? 



Ma 



1 + 



2Ml(M|cos2 2/5 + esin2 /?) 
Mfsin^ 2/? + ecos2/? 



2Ml 



2Ml 



(3) 



The h and if boson masses are displayed in Fig. 2 as a function of the pseudoscalar mass 
for several values of tan/5 = 1.1, 1.6, 5 and 50 and for mq;=250 GeV and 1 TeV. In the 
case of h, the decoupling limit ~ M™^^ is reached very quickly for large values of 
tan/3 [already for Ma ~ 110 GeV] and the maximal h mass is large, up to M™^^ ~ 130 
GeV. For small tan/? values, the maximum h mass is rather small for tan/3 = 1.1; this is 
due to the fact that cos 2/3 is close to zero and Mh is entirely generated through radiative 
corrections. The approach to the decoupling limit is rather slow, and for tan/3 = 1.6, the 
value M^^^ ~ 80-100 GeV is reached only for Ma — 500 GeV. In the decoupling limit, 
the heavy CP-even Higgs particle becomes degenerate with the pseudoscalar, Mh ~ Ma- 
Similarly to h, this occurs very quickly for high tan (3 and slowly for low tan /3 values. 

The charged Higgs boson mass is not affected by the large radiative correction eq. (1) 
and does not depend on tan /3, it is simply given by 



Mh± 



Ma 



1 + 



Ml 



nl/2 



Ma for Ma > M, 



w 



(4) 



It is shown in Fig. 2 together with the h/H masses. Finally, the mixing angle angle a 
which also receives large radiative corrections 



tan 2a; 



M\ + Ml 
^''''^^Ml -M| + e/cos2/3 ' 



TT 

-- < a < 



(5) 
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reaches the values — 7r/2 in the decouphng hmit. 

The two-photon decay width of a CP-even Higgs particle ^ = h,H can be written as 



r($ ^ 77) 



(6) 



128v^7r3 

where the scaling variable is defined as Xj = M|/4m^ with the mass of the loop 
particle. While in the SM one has only contributions from the W boson and heavy 
fermions, in the MSSM additional contributions are provided by the charged Higgs boson, 
the two chargino states and the scalar partners of the fermions; Fig. 1. Factorizing the 
reduced couplings of these particles to the Higgs boson and to the photons, the amphtudes 
of the various contributions read [H 

= g<i>ww Fi{tw) 
Af = N,Q)g^ffFy,{Tf) 



M? 



Ah± = g<S>H+H- Foi'^H±] 

Xi 



fi 

with Nc the color factor and Qf the electric charge of the (s)fermion in units of the proton 
charge. With the help of the function /(r) defined by 

arcsin^x/r r < 1 



the spin 1, 1/2 and spin amplitudes are given by 

Fi(r) = [2r2 + 3r + 3(2r-l)/(r)]/r2 

Fi/2{r) = -2[r + (r-l)/(r)]/r^ 

Fo(r) = [r-/(r)]/r^ (9) 

The amplitudes are real if the Higgs mass is below the particle threshold, Mq, < 2mj, 
while they are complex above this threshold. In the limit of heavy loop masses, r -C 1, 
these amplitudes reach the asymptotic values 

Fi ^ +7 , Fi/2 ^ and Fq ^ (10) 

Note that while the W and fermion loops give finite contributions in the asymptotic limit, 
the contributions of the charged Higgs boson, the charginos and the sfermions vanish in 
the large loop mass limit since the amplitudes Ai are damped by the heavy masses. 
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3. Higgs production at 77 Colliders 

The two-photon width of Higgs bosons can be directly measured at 77 colhders, with 
the photons generated by Compton back-scattering of laser beams from electron beams 
The electron and laser beam polarizations can be chosen such as to tune the photon 
energy spectrum and produce a peak at a fixed energy. The energy of the 77 collider can 
be as much as ~ 80% of that of the original e^e~ collider [|1^, In the following, we will 



discuss briefly the production of the SM Higgs boson in the intermediate mass range at 
photon colliders; we will consider the various backgrounds and estimate the precision with 
which the 77 width could be measured. For this purpose we will follow Ref. |]19 

from which the collider configuration, like beam polarization, conversion distance, etc., 
is taken from. We will assume that the Higgs boson mass Mh is already known and 
therefore, choose the option in which the beam energy is tuned for the 77 luminosity to 
peak at Mh- We also adopt the option where the handedness of the electron/positron 
beams and laser photons are opposite in order to enhance the = partial wave in 
which the Higgs boson signal occurs. 

For masses below < 130 GeV, the Higgs boson will dominantly decay into bb pairs; 
decays into charm quarks, r leptons and gluons occur at the level of a few percent and for 
masses close to Mh ~ 130 GeV, the WW decay mode becomes important and reaches a 



branching ratio of ~ 30% |^. The Higgs boson is extremely narrow, with a total decay 
width below ~ 10 MeV. Being induced by loops, the 77 width is very small, of 

order of a few 10 keV in the mass range 80 GeV < Mh ^ 130 GeV. For a complete 



discussion, see Ref. from which we take the inputs for masses and couplings. 



Since the process we are interested in proceeds through a very narrow resonance, the 
detector accuracy when comparing the 77 —>■ — > bb signal and the background should 
be taken into account. A simple way to obtain the effective signal and backgrounds 
consists of introducing a gaussian smearing of the two-photon invariant mass W, 

da"^ rymV^-+:^ , 1 f (W' - WY 1 dL 



m ' ' Jmx V2^s [ 252 J dw' 

and selecting events within a bin of invariant masses M/^± A. In the previous expression, 
Leff and ymy/Se+e- are the effective luminosity and the maximum energy of the 77 collider; 
S is one sigma of the detector resolution for W. The cross section for the signal process 
77 —>■ bb can be written as 



2r(if0 ^ 77)BR(ifO ^ bb) 

' mT 



asciW) = An'^ ^(1 + X,X2)SiW - Mh) , (12) 



where the helicities of the scattered photons must be such that A1A2 = 1. Inserting the 
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cross section in eq. (11), and selecting the events in the bin Mh ± A, one obtains 



^cff (rf^^Mu) = RiA/5) 



dL 



Jz=0 



Stt- 



W=Mh 



Ml 



(13) 



with R{A/6) being the Gaussian error function, describing the fraction of signal events 
contained in the bin Mfj ± A [for instance, for A = 1.255 one has R ~ 0.75]. 

In the intermediate mass range, 80 GeV < Mh < 130 GeV, the main source of 
background is the continuum production of b- and c-quark pairs, including gluon radi- 
ation which leads to fake two-jet events [^|. The contribution of resolved photons to 
heavy quark production is rather controversial since the partonic distribution functions 



for polarised photons are not yet available; see the discussions given in [|1^, Another 
potentially large background for Mh ~ Mz, is the process 67 — (e)Z — bb which comes 
from the residual electrons that were left over from the Compton scattering; the scattered 
electron is emitted backwards down the beampipe. It can be reduced by removing the 
residual electrons from the interaction region with a strong magnetic field, and this re- 
quires a non zero conversion distance ||23|. Finally, we have 77 Z{ff) bb |l^, ^ 
which would constitute a serious problem that may be overcome using a very peaked 77 
luminosity distribution at ~ Mz to greatly increase the signal-to-background ratio. 

Here, we will only include the continuum qq and qq{g) backgrounds that we calculated 
using the package HELAS . As in Ref . pO| , we have used the following set of experimen- 
tal cuts for the signal and backgrounds: (i) both jets from the q and q should be visible 
in the detector: | cos 6*5^5! < 0.7; (ii) the gluon jet should escape detection: | cos 6*^1 > 0.9; 
(iii) the jets should be clearly isolated: mfj/se+e- > 0.02; (iv) the missing pt and the 
aplanarity due to missing gluon should be small: < 10 GeV and + — 7r| < 0.02. 
For the detector accuracy, we also employ the same resolution as in Ref. ||2^ : 6 = A GeV 
and A = 5 GeV for half of the width of the selection interval. The effective cross sections 
for a tuned energy O.S^Sg+g- = Mh are given in Table 1 for three choices of the Higgs 
boson mass Mh = 80, 105 and 130 GeV. 



Process 


Mh = 80 GeV 


Mh = 105 GeV 


Mh = 130 GeV 


77 ^ bb 


67.9 


73.0 


62.6 


77 bb 


18.1 


7.13 


3.40 


77 — > cc 


240 


99.4 


49.2 


77 bb{g) 


0.13 


0.05 


0.02 


77 ^ cc(^) 


1.72 


0.72 


0.36 



Table 1: Effective cross sections [in fb] for the signal and the backgrounds for a tuned energy 
0.8^Se+e- = Mh and the luminosity distribution of Fig. 16 of Ref. [lI 



For the bb final state, we will assume a detection efficiency of 50% [which should be 
achieved in the future by micro-vertex detectors] with a 5% contamination from cc final 
states. One the other hand, only the events where both quarks decay hadronically should 
be collected to estimate the 77 invariant mass; the hadronic decay branching ratios of b- 
and c-fiavored hadrons are 75% and 82% respectively. Multiplying the tagging efficiencies 
times the square of the hadronic branching ratios, we obtain the corrected effective cross 
sections in Table 2. As can be seen, the radiative background is completely negligible and 
the signal cross sections are much larger than the backgrounds, especially for high Higgs 
boson masses, leading to a large statistical significance for the Higgs boson signal. 



Process 


Mh = 80 GeV 


Mh = 105 GeV 


Mh = 130 GeV 


77 ^ //o ^ bb 


19.1 


20.5 


17.6 


77 —>■ bb 


5.09 


2.00 


0.96 


77 — i> cc 


8.07 


3.34 


1.65 


77 bb{g) 


0.037 


0.014 


0.006 


77 cc{g) 


0.058 


0.024 


0.012 


Total Background 


13.3 


5.38 


2.63 


Signal /Background 


1.44 


3.81 


6.69 


Stat. Significance 


16.6 


27.9 


34.3 


Sensitivity to F^^ 


9.4% 


7.8 % 


8.1 % 



Table 2: Corrected effective cross sections [in fb] to include tagging efficiency and the invariant 
mass reconstruction. An integrated luminosity of 10 fb~^ is assumed for the evaluation of the 
statistical significance and the sensitivity to the two-photon width of the Higgs boson. 

The measurement of F(if° — 77) x BR(if° 77) will follow from eq. (|13]) if the 
luminosity, the tagging and mass reconstruction efficiencies as well as the Higgs boson 
mass are precisely known. Assuming that BIi{H^ 77) is given by the SM and that 
the uncertainties in all the previous quantities are negligible, the statistical error in the 
r{H^ —>■ '-f'-f) determination is 

AT _ 1 y/STB 

r vieff s 

which is about 10% for an effective luminosity L^s = 10 fb~^. Increasing the luminosity 
and improving the 6-tagging efficiency and purity as well as the reconstruction of the bb 
invariant mass would enhance the sensitivity to the /177 coupling. 
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4. Loop contributions in the MSSM 
4.1 W boson loop 

Compared to the SM case where gnoww = 1; the W boson amphtude for the hghtest 
MSSM Higgs particle h is suppressed by a factor Qhww = sin(/3 — a). However, in the 
decouphng regime M| 3> M|, the hWW couphng approaches quickly the SM coupling 



ghww = sm{P-a) 



1 o ^Mj. 
8 ^ Mi 



1 - 



2M| cos 2/3 



(14) 



The W boson form factor Aw is shown in Fig. 3a as a function of both in the SM 
and in the MSSM. The MSSM contribution has been obtained by fixing the pseudoscalar 
mass to Ma — 250 GeV and varying the value of tan/3 from tan/3 = 1.1 to 50. The 
difference between the SM and the MSSM contributions is very small, even for low tan /3 
values where the decoupling limit is not completely reached yet for Ma = 250 GeV. This 
is due to the fact that sin(/3 — a) approaches unity very quickly, the difference being of 
0{Ml/Mi). 



4.2 Fermion loops 

Since the $// couplings are proportional to the fermion mass, the contribution of the light 
fermions to the $77 amplitude is negligible. Only the top quark, and to a smaller extent 
the charm and bottom quark, as well as the r lepton, will effectively contribute. Compared 
to the SM case where gm/f = 1, the huu/hdd couplings are suppressed/enhanced by the 
factors 



9huu 
ghdd — 



COS a 
sin /3 
sin a 

cos /3 



1 + - 



1 - 



IMl 
2 Ml 
I Ml 



2 Ml 



cot /3 sin 4/3 



tan /3 sin 4/3 



1 - 



2M|cos2/3_ 
e 

2M|cos2/5 



(15) 



The fermionic amplitudes At and Ajcr ai'c shown in Fig. 3b as a function of M/j, with 
Ma again fixed to 250 GeV. In the SM, the dominant fermionic contribution At is almost 
constant and can be approximated by At ~ NcQj ^1/2(0) = —16/9. It is smaller than the 
W boson contribution and the two amplitudes interfere destructively. In the MSSM, the 
variation with M^ is rather pronounced. This is due to the variation of the coupling gf^u 
since the decouphng hmit is not reached yet for Ma — 250 GeV and small tan /3 values: 
contrary to ghww, the coupling ghuu approaches the decoupling limit slowly, ghuu 
l-0{Ml/Ml). 

For the bottom quark loop, the amplitude Ai, has both real and imaginary parts since 
Mh > 2mij. The real part of A^, calculated with a running mass m6(M|) ~ 3 GeV, 
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is much smaller compared to At as expected^, but is of the same order as A^. since the 
latter is not penalized by the charge factor. The imaginary parts Im(A{, t-) are larger than 
Re(y4f,T-), but since they do not interfere with the dominant A-\^ and At contributions, 
their effect on the $ — > 77 width is rather small. Note that the difference between the 
SM and MSSM is still rather large for Ma = 250 GeV, however, this difference will hardly 
be noticed in r(/i — > 77) since the contributions of the b and r loops are small. 

Finally, we note that the QCD corrections to the dominant top quark loop are well 
under control and can be included by simply multiplying the Born amplitude by a factor 
(1 — as/vr). The QCD corrections to the b quark loop do not exceed the level of a few 
times ag/n if the running b quark mass at a scale M$/2 is used in the Born amphtude; 



for more details on these corrections, see Ref. |12| 



4.3 Charged Higgs boson loops 

In the coupling of the lightest CP-even Higgs particle h to charged Higgs bosons, large 
radiative corrections which cannot be mapped into the mixing angle a will appear. Re- 
taining again only the leading correction, the Qhh+h- coupling is given by |2^ 



. ,^ . cos2/5sin(/? + a) e cosacos^/3 , , 

g,n.H. = M(3-a) + ^ + 2c^^M| sin/3 ^'^^ 



with = 1 — = sin^ 9w In the decoupling limit, the coupling reduces up to (9(e) 
terms, to ghH+H+ 1 - cos^ 2/?/ (2c^). 

The form factor Ah± is shown in Fig. 4a as a function of Mh± for tan j3 = 1.6, 5 and 
50. Because the contribution is damped by a factor 1/M^± for large masses, and also 
because the spin-zero amplitude Fq is small, the charged Higgs contribution to the h'j'j 
coupling is very small. For low masses, Mh± ~ 100 GeV, Ah± can reach values close to 
~ —0.1, but for Mh± > 250 GeV the contribution of the if^ loop is already only a few 
per mille of that of the dominant W boson loop, and is therefore completely negligible^. 



4.4 Scalar lepton and quark loops 

The left- and right-handed scalar partners of each SM charged fermion, fi and f^, mix 
to give the mass eigenstates /i and f2- The mixing angle is proportional to the fermion 
mass and is therefore important only in the case of the top squarks for the scalar 

^In the MSSM, however, far from the decouphng hmit and for large values of tan (3, the amplitude Af, 
can be very large since the coupling g/jf,;, tan/3 is strongly enhanced. 

•^Note that in two-Higgs doublet models, charged Higgs boson loops will provide the only additional 
contribution to the /177 coupling. Since this contribution is very small, discriminating between this model 
and the SM in the decoupling regime using the ^,77 coupling will not be possible. 
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partners of light fermions the current eigenstates are identical to the mass eigenstates^. 
In this subsection, we will discuss the contribution of slepton and the scalar partners of 
the light quarks only, the contributions of top squark loops will be discussed separately 
later. 

The reduced couplings of the h boson to the left-handed and right-handed partners 
of light fermions, are given by 

9hfUL = (4 - <5/s^) sin(/3 + a) 

dhfnfn = QfSwsm{p + a) (17) 

with // = ±1/2 and Qj the weak isospin and the electric charge of the fermion /. In the 
decoupling limit, one has sin(/? + a) ^ — cos 2/9. 

The contributions of the slepton and squark [except for top squark] loops are shown 
in Fig. 4b as functions of the masses and for the three values tan/3 = 1.6, 5 and 50 
with Ma fixed to 250 GeV. We have summed over all slepton and squark [except stop] 
contributions, and used common masses mj and rriq, which is approximately the case in 
SUSY-GUT models. The form factor Aj is approximately equal to Ah± except that the 
trend for various tan/3 values is reversed. The contribution of the squark loops Aq has 
almost the same magnitude as the contribution Ai, but is of opposite sign. As in the case 
of the charged Higgs boson, slepton and squark loop contributions to the /i — 77 decay 
width are very small: in the decoupling limit and for loop masses above 250 GeV, they 
do not exceed a few per mille and can be safely neglected. 

4.5 Top squark loops 

Due to the large value of the top quark mass, the mixing between the left- and right- 
handed scalar partners of the top quark, ii and in, can be very large. The mass eigenstates 
ti and ^2 are obtained by diagonalizing the mass matrix 

_ ( + COS 2/3(1 - |s^) M| mtm^^ 

* Y mtm^^ mf^ + + | cos 2/3 M| 

where the left- and right-handed scalar masses m^^ and m^^ are generally assumed to 
be approximately equal to the common mass of the scalar partners of light quarks rriq. 
In terms of the soft-SUSY breaking trilinear couplings At and the Higgs-higgsino mass 
parameter /i, the off-diagonal term m^^ reads 

m^^ = At - fi cot f3 (19) 

^The mixing in the sbottom sector can also be sizeable for large values of tan/3. We have checked 
explicitely that this mixing will not affect significantly the numerical results compared to the no-mixing 
case, if the value of the off-diagonal entry in the sbottom mass matrix is not prohibitively large. 
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The top squark masses and the mixing angle are then given by 



1 r 



r — m. + - 
iM 2 



2 I 2 



(20) 



sin 26, 



2mtm\ 



LR 
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COS 26 f 



ti 



2 

mf 



(21) 



The couphngs of the h boson to top squarks in the presence of mixing are given by 



Qhhh ~ 



- sm{a + P) 

rrit sin 26t 
^ 2M| 

1 sin(a + P) 

2 2 
rut sin 26 1 

2M| 



cos 6lt - -s^yCos26't 



cos a sma 

——^^t + — -^/^ 
sm p sm p 



4 

sin^^t + -s^cos2^t 



cosa , sma 
sm (i sm p 



cos a m 



sin (3 M^z 



cos a m| 
sin /? M| 



(22) 



In the decouphng hmit, these vertices reduce to 



9hhh = -^cos2/? 



^fet2t2 = -2 COS 2/3 



cos 6't - -s^cos26't 



4 

sin^ 6't + -s^ cos 26't 



Ml 

2 

Mf 



+ -sin2^t 



2 M| 



(23) 



TOg, the only parameters which enter the contribu- 
77 decay width in the decoupling limit are [that we 



Assuming as usual that m^^ 
tion of the i loops to the h 
will trade against m^J and m^^. There is also a dependence on tan/? which arises from 
the Higgs coupling to top squarks and from the mixing angle since the stop mass matrix 
contains also a small cos 2/3 term. However, this dependence on tan/3 is rather small. 

In Fig. 5a, we show contour plots in the (mf^, m^-J plane for which the contribution Ai 
[which includes the amplitudes of both top squarks] to the h^i^ coupling is \Ai\ = 2, 1, 0.5 
and 0.2. For the sake of convenience, we also display in Fig. 5b, contours in the (m^^, m^J 
plane for fixed masses of the heavy top squark and the scalar partners of the light squarks. 
To have a better insight on the various contributions, we show in Fig. 6 three dimensional 
plots of the two separate top squark amplitudes A^^ and A^^. 

The amplitude A^ is symmetric for positive and negative ^ values because for large 
mf^, the hit coupling is dominated by the sm26t'mf^ term and sin2^t is proportional 
to mf^; for small mf^ the dominant piece of the hit couphng is proportional to mj. To 
discuss the effect of the mixing, it is convenient to divide the parameter space into three 
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regions: intermediate |m^^| values around the region delimited by the contour ^4^ = 0, 
large and small |m^-^| values away from this contour. 



For large \m^^\, the contributions are large and positive; for light enough top squarks, 
~ 100 GeV, they can reach the value ~ 2 for \m^^\ ~ 1 TeV, therefore almost 
canceling the top quark loop contribution At. For a given m^^, is larger for higher 
values of m^^ because in this case, the coupling htt ~ m^^ is strongly enhanced. For 
large m^^, the two top squarks will have comparable masses [see Fig. 5b], and since the 
signs in the dominant component of the htiti and /it2^2 couplings are opposite, the two 
amplitudes will partly cancel each other; Fig. 6. 

For small there is a region [the small "menhir" around m^^ = 0] where no 

solution for m^^ < rrit is allowed when diagonalizing the mass matrix; however, this 
region is already excluded by CDF/DO data from the negative search of scalar partners 



of light quarks with masses rrig < 150 GeV [^; Fig. 5b. The amphtudes in this region 
are negative since the dominant component of the htt is now proportional to and 
has opposite sign compared to the dominant off diagonal coupling when m^^ is large. A^ 
decreases with increasing top squark mass as expected, and can reach the almost maximal 
value At ~ —0.5 for m^^ < 250 GeV; however, most of this region is again ruled out by 
the CDF/DO bound rrig > 150 GeV as shown in Fig. 5b. 

For intermediate values of \m^^\, there is a balance between the two components of 
the htiti coupling which tend to cancel each other, and the two contributions Af^ and 
Af^ which become comparable [since for relatively large m^-^, ti and ^2 have comparable 
masses. Fig. 5b, and A^^ is small] and interfere destructively. At some stage, the two 
contributions cancel each other leading to the contour A^ = of Fig. 5a. As one can 
see, for top squarks not much heavier than ~ 250 GeV, one can have contributions of the 
order of 10% or more to the /i — > 77 decay width if the off-diagonal entry in the stop 
mass matrix is large, Imf"^! > 1 TeV. 

Since top squark contributions can be relatively large, we have derived the two-loop 
QCD corrections to the scalar quark loops in the limit of heavy squarks, using low energy 
theorems. The derivation of the result is done in the Appendix. The effect of the QCD 
corrections is to shift the value of the Born form-factor Ag by an amount 



A _ _ /iBorn 



8 as 

3 71 



(24) 



The correction is about three times larger than in the case of quark loops, and has opposite 
sign. It is of the order of ~ 10%, and therefore well under control. 
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4.6 Chargino loops 



The masses of the two chargino states depend on tan /3, the gaugino mass parameter M2 
and the Higgs-higgsino mass parameter 11: 



1 



m 



XI, 2 



Ml + /i' + 2M^ 



Ta/ (M| - /i2)2 + 4M^ cos2 2/5 + 4M2,(M| + ^j? + 2M2H sin 2/5 (25) 



The chargino couphngs to the hghtest Higgs boson h are given by [pT 



^?»X^Xi 



COS a cos sin 9 + sin a sin 9^ cos 6* 



X2 X2 



COS a cos 9 sin 6'"^ + sin a sin ^ cos 



(26) 



with e = sign(M2/i — sin 2/3) and the angles 6*^ according to 



tan 26'" 
tan 29^ 



2y/2MwiM2 cos /? + /i sin /?) 



Ml 



2M2, cos 2/5 



L2 /i ^Ji'-'vi^' 

2v/2Mh/(M2 sin/? + /i cos /?) 



Ml 



/i2 + 2M^ cos 2/5 



(27) 



Note that in the decoupling limit, one has cos a = sin/5 and sin a = — cos/3. 



The contribution of the chargino loops to the /177 coupling is shown in Fig. 7 in 
the (M2,/i) plane for the two values tan/5 = 1.6 and 50. Contours for \A^\ = 0.5,0.3,0.2 
and 0.1 as well as the region of the parameter space for which the lightest chargino mass 
is larger than 90 GeV [a value below which the charginos will be found at LEP2] and 250 
GeV [which will be probed at a 500 GeV e^e~ collider], have been drawn. 

The chargino contributions to the h'-ff coupling are much larger than those of charged 
Higgs bosons, sleptons and the scalar partners of light quarks. This is due to the fact 
that for heavy particles, the amphtude F1/2 —4/3 is larger than Fq —1/3, and also 
because the chargino contribution scales like ~ V'^^x contrary to the amplitudes for 
scalars which scale like Ai ~ 1/mf. For high tan/5 the contributions are positive, while 
for low tan/5 the amplitudes follow the sign of /z. 

The largest contributions are obtained for small values of tan/5. For chargino masses 
very close to the LEP2 limit, m-^ ~ 100 GeV, one can have to > 0.5 inducing contri- 
butions to r(/i — > 77) which exceed the 10% level. For chargino masses above 250 GeV, 
the maximum A^ contribution will be below ~ 0.2 for both values of tan/5, altering the 
total h ■—' ■j'j width by less than 10% percent. 
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5. Conclusions 



We have analyzed the contribution of charged Higgs bosons, charginos, sleptons and 
squark loops to the coupling of the lightest neutral Higgs boson to two photons in super- 
symmetric theories. Our aim was to determine the region of the MSSM parameter space 
in which one is still sensitive to the additional SUSY loops, although no SUSY particle 
has been produced directly at the LHC or at a future e~^e~ collider with a cm. energy of 
y/s = 500 GeV. We focussed on the decoupling limit where all the additional MSSM Higgs 
bosons are very heavy. In this limit, the h boson will have practically the same properties 
as the standard Higgs particle, and the two-photon decay could be used to discriminate 
between the SM and the MSSM Higgs sectors. Our conclusions are as follows: 

The contributions of charged Higgs bosons, sleptons and the scalar partners of the 
light quarks including the bottom squarks are extremely small. This is due to the fact 
that these particles do not couple to the Higgs boson proportionally to the mass, and the 
amplitude is damped by inverse powers of the heavy mass squared; in addition, the scalar 
loop amplitude is much smaller than the dominant W amplitude. For masses above 250 
GeV, the effect of scalar particles [with the exception of the top squark] on the /i — > 77 
width does not exceed one precent level and can therefore be neglected. 

The contribution of the charginos to the two-photon decay width can exceed the 10% 
level for masses close to ~ 100 GeV, but it becomes smaller with higher masses. The 
deviation of the r{h ^ 77) width from the SM value induced by charginos with masses 
mv = 250 and 400 GeV is shown in Fig. 8 function of M2 [/i is fixed by m^] for 

tan/3 = 1.6 and 50. For chargino masses above > 250 GeV [i.e. slightly above the 
limit where charginos can be produced at a 500 GeV e^e~ collider], the deviation is less 
than ~ 8% for the entire SUSY parameter spaceQ. The deviation drops by a factor of two 
if the chargino mass is increased to 400 GeV. 

Because its coupling to the lightest Higgs boson can be strongly enhanced, the top 
squark can generate sizeable contributions to the two-photon decay width. For stop 
masses in the ~ 100 GeV range, the contribution could reach the level of the dominant 
W boson contribution and the interference is constructive increasing drastically the decay 
width. For ti masses around 250 GeV, the deviation of the h ^ 'jj decay width from the 
SM value can be still at the level of 10% for a very large off-diagonal entry in the stop 
mass matrix, m^^ > 1 TeV [Fig. 8b]. For larger masses, the deviation drops ~ l/'^f^ 
and the effect on the decay width is below 2% for m^^ ~ 400 GeV even at mf'^ ~ 1 TeV. 

^The maximum deviation is obtained for M2 values slightly above the lightest chargino mass. The 
reason is that the Higgs boson prefers to couple to mixtures of gauginos and higgsinos, and in this region 
the hxiXi coupling is maximal. For larger M2 values, xt ^ pure higgsino while for smaller M2 values 
it is a pure gaugino and the Higgs coupling is therefore small. 
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For small values of m^^, the deviation does not exceed —8% even for a light top squark 
r/ij^ ~ 250 GeV. For this stop mass value, we have cut out the region m^^ < 200 GeV 
since there, the scalar partners of light quarks will have masses smaller than rriq ~ 250 
GeV. 

Thus, the only way to have a contribution to the two-photon decay width of the h 
boson that is larger than 10% is from a rather light stop squark m^^ < 300 GeV with 
extremely large couplings to the Higgs boson, m^^ > 1 TeV. However, for m^^ ~ 300 
GeV and m^^ < 1 TeV, the scalar partners of the light quarks will have masses around 
rriq ~ 500 GeV and therefore should be observed at the LHC. The only way to have a 
light top squark, m^^ ~ 300 GeV, while the other squarks are heavier than 1 TeV and 
escape detection at the LHC is to increase m^^ to ~ 5 TeV. For such large mf^ values, 
the trilinear scalar interaction become extremely strong and could lead to dynamically 
favoured minima of the scalar potential where charge is not conserved [0. [A necessary, 
though not sufficient condition to avoid these false vacua is to choose m^^ < 3mg; see 
Ref. [0].] Furthermore, the hit coupling, (yf^^^- ~ mtm|j:j/(2M|), becomes very large and 
perturbation theory is endangered. 

In summary: only chargino and top squark loops can lead to a sizable difference 
between the two-photon decay width of the lightest SUSY and the standard Higgs bosons 
in the decoupling limit. Charginos with masses above the production threshold of a 500 
GeV e"'"e~ collider will induce contributions which are smaller than a few percent in the 
entire SUSY parameter space. Top squarks can induce contributions which exceed the 
10% level only if they are just slightly heavier than ~ 250 GeV and if the off-diagonal 
entry in the stop mass matrix is very large, mf'^ > 1 TeV. In the region of parameter 
space where both charginos and top squarks are hght, the sum of the two contributions 
can exceed the 10% level. 

One therefore needs a measurement of the Higgs coupling to two photons at 77 colliders 
with an accuracy better than 10% to discriminate between the standard and the minimal 
SUSY Higgs scenarios if this discrimination could not be achieved in the e~^e~ option of 
the collider with a cm. energy in the 500 GeV range, or if scalar quarks have not been 
observed at the LHC. 
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Appendix: QCD corrections to scalar quark loops 

The calculation of the QCD corrections to scalar quark loop contributions to the Higgs- 
two photon coupling can be calculated by extending low-energy theorems [l^ to scalars at 



the two-loop level [|12|, ^ . For a CP-even Higgs boson H with a mass, Mh -C 2mg which 
is the case here, these theorems relate the matrix elements of the squark contributions to 
the if77 vertex to the photon two-point function. In this Appendix, we will use these 
theorems to derive the QCD corrections to the scalar loops due to pure gluon exchange; 
this part of the correction is expected to be the dominant one. 

Denoting the matrix element of one squark contribution to the photon self-energy by 
A^g(77), and the corresponding matrix element with an additional light Higgs boson by 
A^q(77/7), one has at lowest order 

MV^ = {V2G,f" m^dM^/dniQ (A.l) 

with Cq and the squark charge and its reduced coupling to the Higgs boson. To 
extend this relation to higher orders, one has to replace all quantities by their bare values, 
differentiate with respect to uiq and then perform the renormalization. In the case of 
pure gluon exchange, the differentiation with respect to the bare squark mass can 
be rewritten in terms of the renormalized mass itlq. A finite contribution to the QCD 
corrections arises from the anomalous mass dimension 7g 

m%-^ = ^ -— (A.2) 

^ om^ 1 + 'Jq orriQ 

The remaining differentiation with respect to the renormalized squark mass of the photon 
two-point function leads to the squark contribution to the /3 function, (3q. The final result 
for the squark contribution to the Higgs-two-photon coupling can be expressed in terms 
of the effective Lagrangian 

Ci = {V2G,f\lf^F>^^F,.H (A.3) 

The QCD corrections to the squark loop are then fully determined by the anomalous mass 
dimension of the squarks 7q = 4a<j/ (Svr) and by the squark contribution to the /3 function 
/Sq/a = 2a/'7r(l + Aas/ir) ||3^. This results into a final rescaling of the lowest-order 
Lagrangian by factor 



3 vr 



(A.4) 



This has to be compared with the case of heavy quark loops where the QCD correction 
gives a rescaling factor 



vr 



(A.5) 
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Figure Captions 

Fig. 1: Feynman diagrams contributing to the Higgs decay into two photons in the MSSM. 

Fig. 2: The masses of the three MSSM Higgs bosons /i, H and as a function of the 
pseudoscalar mass Ma for tan/j = 1.1, 1.6, 5 and 50; the common squark mass is 
set to ruq = 250 GeV (a) and 1 TeV (b) and the squark mixing is neglected. 

Fig. 3: The amplitudes for the contribution of the W boson loop (a) and of the t, b loops 

(b) as a function of Mh in the SM [dashed lines] and MSSM [sohd lines]; Ma is fixed 
to 250 GeV. We have used rriq — 250 GeV and neglected squark mixing. 

Fig. 4: The amphtudes for the contribution of the charged Higgs boson loop (a) and of 
the slepton and squark (except stop) loops (b) as functions of the loop masses for 
tan/3 = 1.6,5 and 50. We have neglected squark mixing and for H"^ and slepton 
loops we used rriq — 250 GeV; for sfermion loops we have set Ma — 250 GeV. 

Fig. 5: Contours in the {m^^,mij plane, for which: (a) the contribution of the stop loops 
to the /i77 coupling is \Af \ — 2, 1, 0.5 and 0.2, and (b) the heavier top squark mass 
rrif^ and the common mass of the scalar partners of light quarks rriq. 

Fig. 6: The separate contributions of the lightest (a) and the heaviest (b) top squark loops 
to the form factor Af as a function of mf^ and m^^ . 

Fig. 7: Contours in the {M2, n) plane for tan/3 = 1.6 (a) and tan/3 = 50 (b) for which the 
contribution of the chargino loops to the h'yy coupling is \A^\ = 0.1, 0.2, 0.3 and 
0.5. Also included are the contours for which the lightest chargino mass is = 90 
and 250 GeV. 

Fig. 8: The deviations of the SUSY Higgs coupling to two photons from the Standard Model 
value [in %] for two values of tan/3 = 1.6 and 50 and the loops masses rrii = 250 
and 400 GeV. (a) Deviations due to the chargino loops as a function of M2 for both 
signs of 11, and (b) deviations due to the top squark loops as a function of m^^. 
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Fig. 5a 




Fig. 5b 
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Fig. 6 
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